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We introduce in this work a class of relativistic models for 
nuclear matter and neutron stars which exhibits a parame- 
terization, through mathematical constants, of the non-linear 
meson-baryon couplings. For appropriate choices of the pa- 
rameters, it recovers current QHD models found in the lit- 
erature: Walecka and the Zimanyi and Moszkowski models 
(ZM and ZM3). For other choices of the parameters, the 
new models give very interesting and new physical results. We 
have obtained numerical values for the maximum mass and 
radius of a neutron star sequence, red-shift, hyperon popu- 
lations, radial distribution of particles, among other relevant 
static properties of these stellar objects. The phenomenology 
of neutron stars in ZM models is presented and compared to 
the phenomenology obtained in other versions of the Walecka 
model. We have found that the ZM3 model is too soft and 
predicts a very small maximum neutron star mass, ~ 0.72Mq . 
A strong similarity between the results of ZM-like models and 
those with exponential couplings is noted. Sensibility of the 
results to the specific choice of the values for the binding 
energy and saturation density is pointed out. Finally, we dis- 
cuss the very intense scalar condensates found in the interior 
of neutron stars which may lead to negative effective masses. 

26.60.+C, 21.65.+f, 95.30.Cq, 97.60.Jd 



I. INTRODUCTION 

Since the discovery of the first pulsar in 1968 fif] and 
its identification as a rotating neutron star, the structu- 
re, composition, dynamics and evolution of these astro- 
physical objects became important themes of theoretical 
and phenomenological research. According to an early 
suggestion 0, neutron stars evolve from an initially hot 
protoneutron star that forms in the collapse of a massive 
star in the supernova phenomenon. At densities exceed- 
ing that of nuclear matter, important static properties 
of a neutron star as the mass-radius relation, the crust 
extent, the distribution of the stellar moments of inertia 
and the central density, may be determined by its equa- 
tion of state (EOS) §. 

During this period, there has been a continuous en- 
hancement concerning relativistic microscopic calcula- 
tions of the EOS of neutron stars, improving our un- 
derstanding of the structure of these stellar objects. In 
particular, any theory must at least account for a neutron 
star as massive as the most massive observed pulsar; the 



knowledge of pulsar masses provides a very important 
constraint on the theory. 

More recent calculations based on relativistic proper- 
ties of nuclear matter at high densities indicate that the 
equations of state are considerable stiffer than those pre- 
dicted by non-relativistic approaches. As a result, the 
mass of a neutron star is believed to be at least as large 
as 1.6- 2.1M Q . 

From the theoretical point of view, quantum chromo- 
dynamics (QCD) represents the most profound descrip- 
tion of the strong interaction and would be the ideal tool 
for neutron star applications. However, the highly non- 
linear behavior of QCD at the hadronic energy scales in- 
hibits any practical calculations leading most theorists to 
search for phenomenological descriptions of the structure 
of nuclear matter. One of these alternative approachs is 
quantum hadrodynamics (QHD) [||, a relativistic quan- 
tum field theory based on a local Lagrangian density 
which uses baryon and meson fields as the relevant de- 
grees of freedom. This model provides consistent the- 
oretical framework for describing such a relativistic in- 
teracting many-body system and, based on it, Glenden- 
ning presented Q a very comprehensive treatment of the 
matter in neutron stars using an extended version which 
included leptons and the fundamental baryonic octet. 

Alternative versions of the Walecka model, namely the 
Boguta-Bodmer (BB) || and the Zimanyi-Moszkowski 
(ZM) models |7]||, were developed to improve the de- 
scription of the nucleon effective mass, M*, (too low) 
and compression modulus of nuclear matter, K, (too 
high) as attained with the original approach. Boguta 
and Bodmer introduced cubic and quartic scalar self- 
interactions in the Lagrangian while the ZM models have 
heightened the delineation of these quantities by replac- 
ing the Yukawa scalar coupling term by a derivative cou- 
pling contribution. This derivative coupling may be in- 
terpreted alternatively as a phenomenological coupling 
between the scalar neutral meson and the nucleon fields 
through the introduction of a baryon density dependence 
in the scalar, vector and isovector coupling constants of 
the theory |uj . 

In this work we analyze the structure of neutron stars 
by introducing a QHD Lagrangian with a parameterized 
meson coupling contribution flllfl . This phenomenological 
approach contains high order self-coupling contributions 
of the meson fields and permits in particular to restore 
the results obtained with Walecka, ZM and ZM3 models 
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by making suitable choices of the values of the mathe- 
matical parameters of the theory. The control on the 
analytical form of the couplings allows us to investigate 
other values of these parameters which give new physical 
results. By extending the formalism to include hyperons 
and leptons, we investigate several static bulk proper- 
ties of neutron stars using the Walecka, ZM and ZM3 
models, the two later being applied to this problem for 
the first time. As we have a class of models we are able 
to relate nuclear matter and neutron star properties. In 
particular, we have found that some of the studied mod- 
els describe very strong scalar fields in the interior of 
the neutron stars, leading to a negative nucleon effective 
mass. 



II. GENERAL CHARACTERISTICS OF 
NEUTRON STARS 

A. Electrical Charge Neutrality 

According to Glendenning in ||, la raison d'etre of 
neutron stars is to be neutral. To prove this assertive, 
one should consider that neutron stars are held together 
by the gravitational attraction and take into account the 
balance between the repulsive Coulomb force acting on 
a charged particle of the same sign as the net charge of 
the star (Z net ) and the gravitational force. Assuming the 
particle to be located at the surface of the star, it will be 
expelled out unless the gravitational force overcomes the 
Coulomb force. For the proton, the corresponding limit 
on the net positive charge of the star is Z ne t/A ~ 10~ 36 
fHjf . For the electron, this limit would be reduced by the 
factor m e /m p . Hence, the net charge per charged particle 
is practically zero. This result leads to the conclusion 
that a neutron star is electrically neutral. 

B. Chemical Equilibrium 

In the evolution of the protoneutron star many differ- 
ent reactions can occur. Electric charge and baryon num- 
ber are conserved on a long time-scale in comparison to 
the lifetime of the star. In the core of the protoneutron 
star, the Fermi energy of the nucleons exceed the hyperon 
masses and these particles can be produced in strong in- 
teraction processes with conservation of strangeness in 
reactions such that 

n + n— >n + k + K. (1) 

However, strangeness is not conserved in the time-scale 
of the star since there occurs diffusion of neutrinos and 
photons to the surface of the star and processes like 

K°^2 l; K~ -» fx- + V M , ... (2) 

cannot be reversed anymore and a net strangeness ap- 
pears. 



In this evolution process, the star reaches the chemical 
equilibrium, a degenerate state where, from the point of 
view of its hadronic and leptonic composition, further 
reactions are not possible. As an example, in an ideal 
degenerate system of protons, neutrons and electrons at 
chemical equilibrium, particle levels are filled in such a 
way that neutron beta decay or proton inverse beta decay 
are not energetically favored. 

In general, if one takes into account the fundamental 
baryon octet and lepton degrees of freedom (electrons 
and muons), the following chemical equilibrium equation 
then holds (see Appendix |X|) 

Mi = <?MM« - qe,iHe , (3) 

where i = p,n, A, £ , S + , S~, S°, e, p,; qb,i represents 
the baryon number and q e ^ the electric charge of species 
i. In this way, the conditions for (3 equilibrium can be 
summarized as 

MS = ME = MA = Mn ; 

MS" = MH~ = Mn + Me ; (4) 
M£+ = Mp = Mn _ Me • 

III. BOGUTA-BODMER MODEL 
A. Theory 

In this section we study baryon and lepton populations 
in neutron stars by using the Boguta-Bodmer model (BB) 
with hyperon degrees of freedom. In spite of this study 
having been already done (H|, we reproduce its main re- 
sults as a guide for the development in the next section 
of our new class of non-linear relativistic models. 

The BB model describes the complex fermionic com- 
position of neutron stars as a generalization of the a, lo 
and g theory 

£=^-05 [il^ L - {M B - gaB<j)- 9luBJ^}^B 
B 

-J2 ^BilggBl^r ■ Q^B+Ud^ad^a - m ff 2 er 2 ) 

B Z 

-hM{g a <j) 3 - \c{ga<yf+Y^ - m t )i/>t ■ 

i 

(5) 

This Lagrangian density describes a system of eight 
baryons (B = p, n, A, S°, £+, S°) coupled to 
three mesons (a, uj, g) and two free lepton species {I = 
eT , The scalar and vector coupling constants in the 
theory, , and the coefficients b and c are determined 
to reproduce, at saturation density po = 0.15/m -3 , the 
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binding energy of nuclear matter, B = — X6MeV, the 
compression modulus of nuclear matter, K = 250Af eV, 
and the nucleon effective mass, M* JM = 0.75. In fact, 
the values for these two last quantities are not well estab- 
lished and we have just taken the most used values in the 
literature of the BB model. Additionally, to describe the 
symmetry energy coefficient, — 32.5A/eV r , we deter- 
mine the isovector coupling constant g e . We have found 



9.86 fm 2 



9uj 



5.85 fm 2 



^- ) = 4.80 fm 2 ; b = 0.00103 ; c = 0.0100. (6) 

m, 



In the comparison with the results obtained by ||, 
one should recall that this author has fitted the cou- 
pling constants of the theory with B = — 15.95MeV, 
po = O.U5fmr 3 ,K = 285MeV, M*/M = 0.77 and 
a 4 = 36.8MeV. 

Using the Euler-Lagrange equations, the Dirac equa- 
tion for uniform matter, in momentum representation, 
is 



^ fl (*)=0, (7) 



where M B (a) = Mb — g a B& is the effective mass of the 
baryonic species B. Furthermore, by applying the mean 
field approximation, the u>o, £?03 and a meson field equa- 
tions for uniform static matter are 

E9uiB 
—2?b , 



003 



E 



9qB t 

2 1 3BPB 



m„ 
B e 



+ E 



-bMg G {g a af 



2J B + 1 



2tt 2 



-9aB 



kp - B M B {o) 



Vk 2 



M* B 2 



-.k 2 dk. 



(8) 



In these equations, the baryon source terms have been 
replaced by their ground state values. 

Defining the ratio between meson-hyperon and meson- 
nucleon coupling constants as 

_ 9{cr,uj,g),B 

,g),B 



we have, from equations 
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(10) 
(11) 



E XoB 
2tt 2 



bM(g a a) 2 - c{g a af 
kF - B M* B {a) 



Vk 2 



M* B 2 



k 2 dk] . (12) 



The corresponding equations for baryon number and 
electric charge conservation are: 



and 
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(13) 



(14) 



The baryon chemical potentials, ps (k), correspond to 
eigenvalues of the Dirac equation (0) : 



HB(k) = g^B^o + 9gBQmhB + \J k 2 FB + M B (a) 2 . (15) 

In this expression, I^b is the isospin projection of baryon 
charge states B and kF,B is the Fermi momentum of 
species B. 

The EOS is obtained from the ground state expecta- 
tion value of the time and space components of the diag- 
onal energy-momentum tensor. The energy density and 
pressure of the system are given, in the BB model, by 

£ = ^bM(g a af + -c(g a af + + t^^q 



1 



TO p£03 + E ^2 
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(16) 
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(17) 



We present in the following the results obtained with this 
model. 



B. Results 

In the numerical calculations with the BB model we 
have considered matter with and without hyperons, in 
order to understand how these strange species affect the 
neutron star properties. 



1. Matter with nucleons, hyperons and leptons 

Figure [l], panels a-d, shows baryon and lepton popula- 
tions and field strengths as a function of the total baryon 
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density for two different ratios: \ — -\/2/3 [j and x = 1 
(universal coupling) . From expression ( |15[ ) for the baryon 
chemical potential we can see that the charge term in the 
eigenvalue determines whether a species is charge-favored 
or unfavored; and the isospin term determines whether a 
species is isospin-favored or not. Baryons with the same 
sign of the electric charge as the proton are unfavored; 
baryons with the same sign of its isospin projection are 
favored (notice that g e Qo3 < 0). 

At high densities (p ~ 0.8 /m~ 3 ) the A hyperon be- 
comes the most populous species for the case x — \/2/3 
(panel a). At panel b we see that the electron chemical 
potential reaches a maximum value at ~ 200 MeV and 
begins to decrease due to the reduction of the electron 
population. At p ~ 1.5 fm~ 3 the nucleon effective mass 
is still at 200 MeV. 

In the panel c of the same figure, we see that, to the 
increasing of the meson-hyperon coupling constant (from 
\/2/3 to 1) it corresponds an early emergency of the par- 
ticles. In the comparison with the results of panel a, we 
can see that the leptons have a greater population in this 
case and that the neutron population remains always as 
the most important in the system. The electron chem- 
ical potential and —g g Qo3 saturates around 200 MeV at 
panel d and the nucleon effective mass behaves similarly 
to the case x — >/2/3. 

2. Matter with nucleons and leptons 

In the sequence of the analysis of the results of fig- 
ure [l], panel e shows baryon and lepton populations and 
the corresponding chemical potentials when we exclude 
hyperon degrees of freedom. In this case, since charge 
neutrality is kept only by the p + ,e _ and p~ particles, 
the lepton population increases in the domain of densi- 
ties shown in the figure. Also, we notice the increasing 
of the nucleon effective mass, which means a less intense 
scalar field when compared to the previous cases. This 
indicates that the introduction of hyperons enhance the 
scalar meson condensation. 



3. Neutron star properties 

We are able now to find numerical results for the 
EOS using equations ( |l6| ) and (|l7|). However, this EOS 
corresponds to neutron star matter densities (10 13 — 
10 15 g/cm 3 ) and should be supplemented by EOS's from 
other models for sub- nuclear densities: we adopt the ap- 
proach developed in |l4j] in the density interval 2 x 10 3 — 
lx 10 11 gf/cm 3 and in the range lx 10 11 — 2x 10 13 <?/cm 3 we 



1 This choice is based on the quark counting of the baryons 



use the EOS presented in [[l5|. Combining these EOS's 
with the Tolman-Oppenheimer-Volkoff (TOV) equations 
]l6| , pTt , we may determine the mass of a neutron star as 
a function of its central density. To gain the mass-radius 
relationship, the radius R of a neutron star is obtained 
with the condition that the pressure is null at the sur- 
face of the star, p(R) — 0. We have found values for the 
mass and radius of different neutron star sequences as a 
function of the central density e c . 

To illustrate the importance for neutron star struc- 
ture of the EOS at densities greater than 10 13 g/cm 3 , 
we present in figure |[a the results of the energy density 
as a function of the radial coordinate for three specific 
values of the neutron star mass, for the case without 
hyperons. We can see that the energy density is approx- 
imately constant for heavier stars and is found mostly 
around 10 14 g/cm 3 . 

The conversion of nucleons into hyperons reduces the 
Fermi pressure associated to the baryons, softening the 
equation of state and lowering, as a consequence, the 
maximum mass of a neutron star sequence. In figure ||.b, 
we see the neutron star mass and central density relations 
for the situations analyzed above. The results indicate, as 
expected, that the presence of hyperons causes a diminu- 
tion of the maximum mass of a neutron star sequence. 
Typical results for the mass-radius relation are shown in 
figure pic. 

Figure || shows the behaviour of the radial distribu- 
tion of baryon and lepton populations. The results indi- 
cate that neutron populations are, in general, dominant. 
However, in the case with x — \/2/3 we can see the 
important contribution of the A hyperons in the inner 
regions of the star. Moreover, we have obtained ~ 11km 
for the radius of the star with hyperons and ~ 12km 
without these particles. 

IV. MODELS WITH DERIVATIVE COUPLINGS 

The BB model has two additional coupling constants, b 
and c, associated to self-interactions of the neutral meson 
scalar field. This allows a very good description of two 
important properties of nuclear matter, the compression 
modulus and the nucleon effective mass, which concerns 
the high-density behavior of the equation of state. How- 
ever, some authors argue that the model suffers of a very 
serious problem: the constant c has negative values for 
several entries of M* and K, allowing the energy density 
to become unbounded from below for large values of the 
scalar meson mean field, leading to unphysical behaviour 
& 

In the derivative coupling model, introduced in 1990 
0, the deficiencies of the original Walecka approach 
are eliminated at the cost of making the theory non- 
renormalizable. ZM models have been used in the des- 
cription of static properties of neutron stars |fl9|] , A ex- 
citations in nuclear matter pfj , bulk properties of finite 
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nuclei |T(J, in- medium quark and gluon condensates and 
restoration of chiral symmetr y |^l||, and thermodynamics 
properties of nuclear matter |22jj23| ]. 

The authors of M have presented two additional ver- 
sions of the derivative coupling model. These three mod- 
els are known as ZM, ZM2 and ZM3 §,§. Concern- 
ing the description of static properties of nuclear matter, 
the ZM2 model does not exhibit fundamental differences 
from the ZM model and will not be considered in the 
present study. 

The Lagrangian density of the ZM and ZM3 models 
are: 

C ZM = -tpMijj + m*- 1 ^^^ - g^l^^] 



2 a 2 ) - -u^u^ + -m^uj^, 

(18) 



and 
r 



-ipMijj + rn* 1 i/jij fJi d tJ 'ijj — Qw^la^^ 



-{d^a-m^a 2 ) 



where, 



TO = (1 H ) . 

1 M 



(19) 



(20) 



Expanding to* in expressions ([L8|) and (|19|), in terms of 
the ratio ^ff < 1, we see that the Walecka and the ZM's 
models differ essentially on the replacement of the mini- 
mal coupling Yukawa term g^ipaip, by the derivative cou- 
pling contribution (g CT g/M)-07 M ffl'-0. 

Rescaling in ( |l8| ) and ( p9| ) the nucleon fields in the 
form Ft 



J, 



1 m*i\) , 



(21) 



the Lagrangian densities of the ZM and ZM3 models may 
be written in the general form ^] 



C = $ [i^(d» + ig*u») - (M - g*a)] 1> 

+ ^(d^ad^a - m 2 a 2 ) - -lu^lu'- 11 ' + -m UJ 2 uj f _ l uj f ' 



with: 
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M I 



(1 



M I 



(22) 



(23) 



2 In fact, this rescaling introduces a spurious factor T = 
|(?/>7 (J i/j)c' t '/n(m*(cr)) which does not contribute to the dy- 



namics of the system 
3 Notice that in the ZM models m*M 



Expression (22) reproduces the Lagrangian density of the 
Walecka, ZM and ZM3 models with the following replace- 
ments: 



Walecka : g* -> g a , g* ->• g u ; 

ZM : gl -* m *g<j, gt -> 9u>; 

ZM3 : g* ► m*g a , 9 ^mV 
The resulting field equations are: 

[»7 M 0" - (M - 5 » - <7*7,^1V> = 0; 

cW + to w 2 c^ = 3 ^7^; 

OCT OCT 

which, after the mean field approximation, become 

„2 



- P . 



g<j° 
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for the ZM model and 



g <j da 



* 9 uj 
m —tP, 
mi 



9<J a 



1 dg*a 
da 



tier 



1 dg u 



gag^ 



Oa 



(24) 

(25) 
(26) 

(27) 

(28) 
(29) 

(30) 
. (31) 



for the ZM3 model. A peculiar aspect of the ZM3 model 
is the coupling between the scalar and vector meson 
fields, a kind of coupling which is not present in other 
QHD models (Walecka, BB, ZM). The scalar and vector 
mean field potentials are defined to be 



and 



S = -g>, 



V = g*u • 



(32) 



(33) 



M 



g a a. 



In these models, the expressions for the energy density 
and pressure are equivalent to the corresponding ex- 
pressions ( p"6| ) and ( p"7| ) of the non-linear model, with 
b = c = and the replacements shown in (|24|) . 

The values of the scalar (g a /m a ) 2 and vector (<7 w /to w ) 2 
coupling constants which reproduce the saturation den- 
sity, po = 0.15/to -3 , and the saturation energy, B = 
— 16MeV, in case of the ZM and ZM3 models are shown 
in table || the table also contains the values of the nu- 
cleon effective mass, the compression modulus of nuclear 
matter, the scalar and vector mean fields S and V for 
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both models, and values for the isovector coupling con- 
stant (g g /m e ) 2 which reproduce the symmetry energy 
coefficient 04 = 32.5 MeV. 

From figure |[a we see that the ZM and ZM3 models 
produce a softer equation of state when compared to the 
corresponding results for the Walecka and BB models. 
The ZM and ZM3 models also produce higher values for 
the nucleon effective mass as shown in figure |].b. Fig- 
ures [|a and |5|.b show the behavior of the non-linear cou- 
pling constants (g a * /m a ) 2 and (,g w */m w ) 2 as a function 
of density for the ZM and ZM3 models. Notice that the 
coupling strength decreases with increasing density. 

A interesting comparison among the models considered 
in this section (Walecka, BB, ZM and ZM3) is related to 
their corresponding values for the relativistic coefficient 
defined as: 



R 



P 



k F 



d A k- 



M* 



y/k 2 + M* 



d A k. 



(34) 



According to this definition, for a less relativistic model 
p s becomes closer to p because k <C M. For the models 
discussed in this section we have Rw = 0.931, Rbb = 
0.963, R ZM = 0.970, R ZM3 = 0.959. The results of the 
ZM3 model for the relativistic coefficient are closer to the 
corresponding results of the Walecka model. 

After the work of Zimanyi and Moszkowski, many au- 
thors started to explore extensively the freedom in the 
choice of the meson-baryon interaction. As explained 
below, based on these works we introduce in this paper 
a new class of models which enable us to make direct 
comparisons among the properties of nuclear matter and 
neutron stars. We also employ special cases of this class, 
namely Walecka, ZM and ZM3 models, in the description 
of static properties of these stellar objects. 



V. A CLASS OF NON-LINEAR RELATIVISTIC 
MODELS 

In this section we propose a new class of relativistic 
hadronic models which exhibits a non-linear parameter- 
ization of the intensity of the meson couplings and in- 
corporate some QHD models found in the literature. We 
study, through this comprehensive approach, the influ- 
ence of non-linear meson-nucleon couplings in the nucleon 
effective mass, compression modulus of nuclear matter, 
relativistic coefficient and static properties of neutron 
stars. 



A. Phenomenological Lagrangians 

Koepf et al. p4[ have studied the contribution of the 
term £ ct jv = Mm*(a)'ipij> using the phenomenological 
couplings shown in table |l| The first term in the ta- 
ble corresponds to the Walecka model and the last term 



to the ZM model. Glendenning et al. have analyzed 
a coupling term of the type 



(i + »e£)- 



(35) 



obtaining M*/M = 0.796 and K = 265MeV. 

In these studies the authors have assumed, at first 
order in g a ujM^ a similar expression for the nucleon 
effective mass as introduced by Walecka: M*/M ~ 
1 — g a cr/M. This basic assumption is essential, since 
g a o~ represents the nucleon mass shift in the scalar me- 
son mean field condensate. The experimental results sug- 
gest, at nuclear saturation density, M* /M ~ 0.70, giving 
g a a/M ~ 0.3. Thus, these different models found in the 
literature just add scalar self-coupling corrections terms 
to the corresponding expression of the Walecka model. 

On the basis of the various approaches found in the lit- 
erature, we propose a new phenomenological Lagrangian 
with a non-linear parameterization, through mathemati- 
cal constraints (A,/? and 7 parameters), of the analytical 
form of the meson-baryon couplings: 

B 

B A 



1 



M v -I- -m 2 n . 



where 



9oB — m *\B9v '1 9uB = m ^B9^ 



(36) 



9* e B = m* lB 9g (37) 



and 



= (1 



9<7 0~ 



A,/3,7- 



(38) 



In these expressions, we assume A, (3 and 7 as real and 
positive numbers, since this is the range of best phe- 
nomenology. As discussed in the previous section, es- 
sentially what was done is the introduction of a rescaling 
of the scalar and vector coupling terms of the Walecka 
model. For instance, in case of the scalar contribution 
we have made the replacement 



gao-t/jip -> g^tpo-tp = 



\ L T am; 



-Iplp. 



(39) 



Similar interaction terms may be associated to the vector 
and isovector sector of the Lagrangian density. Notice 
that we have assumed an universal coupling by setting 

9{cr,u,g)B-± 9(a,uj, e )- 



Table [II exhibits the correspondence between this and 
the other models discussed in this work with specific val- 
ues of A, (3 and 7. One of the main intentions of the 
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present study is to consider values of these parameters 
which give better results for nuclear matter and neutron 
star properties when compared to the corresponding re- 
sults of the traditional models discussed in this work. As 
far as we know, the first extensions of the ZM-like models 
to applications to neutron star matter with the inclusion 
of hyperons (A, and leptons was done in p5fl . Here 
we consider these known models as well as intermediate 
values of the parameters of our non-linear coupling to ob- 
tain results for neutron star properties and relate them 
to nuclear matter saturation observables. 

Using the mean field approximation, the field equations 
in our approach become 



[il^ - 5* b7 o^° - (M B - g* B a)]il> B = : 



= (If:) Y, m lBPB\ 



9a 

m„ 



Q, B 
2 



3bpb; 



(40) 
(41) 

(42) 



£(^,,, 



L B 



V 9° 



9a 



E 



G b (<t) 
H B {a) 

g<yg e 



PB 



Hbpb 



where Fb(ct), Gb{<j) and Hb(<j) are given by 



F(a) 



g a <j 

M 



(™*\) 



(a+i)M 



and 



g^guj 



H b (<j) 
g<yg e 



( m *)W+i)//s 
M 



( w *)(7 + D/7 
Mr 



(43) 



(44) 



(45) 



(46) 



From the eigenvalues of the Dirac equation, the Fermi 
energy is: 



MbO) = 9tB^o+9*g B Qo3 I 3B+yk% B +(M B -g* B a) 2 . 

(47) 

The expressions for the scalar and vector potentials (S) 
and (V) are 



and 1, keeping (3 = 7 = 0, we obtain values of S, V, M* 
and K which correspond to the intermediate region of 
values of Walecka and ZM models. Similarly, for values 
of A, (3 and 7 between and 1, we can find intermediate 
results between the Walecka and the ZM3 models. 

Indeed, the range of possible values for the parameters 
of the theory is not very large. Due to the form of the 
general coupling terms (see equation (^8|)), there occurs 
a rapid convergence to an exponential form. Taking A, 
/?, and 7 — ► 00, we have F} 



9o 



g<r ; 

~9q- 



9,. 



(49) 



As we will see later, for A and/or (3 = 7 > 2 the results 
of this model do not strongly differ from the results of 
the model with exponential coupling. 

In this work, we shall consider two cases: 

• Scalar (case S): we consider in this case variations 
of A with j3 — 7 = 0; this case contains the results 
of the Walecka and ZM models. 

• Scalar- Vector (case S-V): we consider in this case 
variations of A, with j3 = 7 = A; Walecka and ZM3 
models belong to this category. 

Notice that the models we are considering may be 
uniquely specified by the A parameter. Walecka model 
belongs to both categories because in this model the 
mathematical parameters A, f3 and 7 are null; but it does 
not present scalar-vector interaction contributions. Tak- 
ing the limit A — > 00 for the S case, we have a model 
with exponential coupling; in the figures we refer to this 
asymptotic model as EXP/S. Similarly, for A,/3,7 — > 00 
we have the asymptotic model EXP/S- V. In the following 
we exploit the nuclear properties of our non-linear class 
of models. 



B. Nuclear Properties 

We determinate the coupling constants g a -,uj,g/'m l7iU j i g 
in this model by following the same procedure presented 
in section III. For each case, S and S-V 
merical values for (/ct/too-, g^/m, 
of A. Thus, we can also determine g e /m g . We also find 
an analytical expression for the compression modulus of 
nuclear matter as a function of the nuclcon effective mass 
(see Appendix |b|) . 

Figure [| exhibits the dependence of the coupling con- 
stants on the A parameter; it is interesting to note 



we obtain nu- 
and M*, as a function 



(48) 



We can see that this model allows some control on the 
intensity of the scalar and vector mesons mean-field po- 
tentials. For instance, to the variations of A between 



4 Notice that this form of the exponential coupling term is 
different from the corresponding one of model 3 shown in table 
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the regular relative behaviour of the coupling constants 
(g a /m a ) 2 and (g^/m^) 2 . The results indicate that the 
scalar model suffers a A dependent saturation process in 
a small range of values for this parameter. On the other 
hand, the results also reveal that the scalar-vector model 
exhibits a wider range of values of A. 

In figure ^ we show results corresponding to the So x 
Vo plane for different values of A. To understand these 
results one should consider expression ( |47j ) for symmetric 
nuclear matter (with k n = k p and the contributions of 
hyperons and leptons turned off) at saturation density: 



Mo = V + ^/k 2 + (M + S ) 



(50) 



Since po = £0/ Po — M + B and using ko/M and So/M 
as expansion parameters, equation ( |50| ) may be approxi- 
mately expressed in the form 



B ~ (V + S ) 



(51) 



This is a very well known and interesting result: the 
Lorentz structure of the interaction leads to a new en- 
ergy scale and the small nuclear binding energy per nu- 
cleon arises from a delicate cancellation between the large 
components of the scalar attraction and vector repulsion 
plus an additional kinetic energy term of a nucleoli with 
non-linear mass M + So (~ M/ (1 — So/M)) at the Fermi 
level. Expression ( pl|) leads to a linear relation between 
So and Vo , explaining the results of figure |7|. 

Results for the compression modulus as a function of 
the A parameter for the S and S-V cases are shown in 
figure @; according to these results, there is a minimum 
of K, in the S case, for A ~ 0.2. In figure ^| we show 
the monotonic behaviour of the nucleon effective mass, 
at saturation density, as a function of A. 

Figure [w| presents the relation between the compres- 
sion modulus and the ratios of the nucleon effective mass, 
K x M*/M, for the S and S-V cases. From the results, 
one can see that to higher values of the nucleon effec- 
tive mass it corresponds lower values of the compression 
modulus: to understand this behaviour one should re- 
member that the repulsive mean-field vector potential 
Vo is proportional to M — M* (see figure [j]) . From 
the results found in the literature p{| , the nucleon effec- 
tive mass and the compression modulus should be in the 
range 0.6 < M*/M < 0.7 and 200 < K < 300 MeV. As 
stressed before, the scalar case leads to a narrow interval 
of values of A while the scalar-vector case has a broader 
range. Accordingly, case S gives reasonable results in the 
range A ~ 0.05 — 0.07 and case S-V better results for 
A ~ 0.16— 0.4. The analysis also reveals the strong sim- 
ilarity between the results of ZM-like models and those 
with exponential couplings. 

The results of figure [Li] show that models with higher 
relativistic coefficients (W and ZM3) have a smaller 
scalar potential S. To understand these results, one 
should consider the definition of R (equation |34|). For 



less relativistic models, R — > 1, since for these models p s 
becomes closer to p because k/M* -C 1. 

We turn now to the high density behaviour of these 
approach by considering its application to the neutron 
star environment. 



C. Neutron Star Properties 

In this section we consider the determination of neu- 
tron stars properties using our class of non-linear models 
with the inclusion of hyperon and lepton degrees of free- 
dom. 



Following the same procedure of section [II, we solve 
a system of transcendental equations taking into account 
chemical equilibrium, baryon number and electric charge 
conservation and the equations for the cr-field. We then 
obtain the EOS for our system. The resulting expressions 
for the energy density and pressure are, again, similar to 
@ and @ but with b = c = and M% = M B - g* B a. 
Combining this EOS with the TOV equations we obtain 
values for static properties of neutron stars (mass, radius, 
baryon composition among others) as functions of the 
central density. 

We present explicitly here the neutron star phenome- 
nology for the Walecka and Zimanyi-Moszkowski models, 
namely the original ZM and also the variant ZM3 model, 
whose results are presented here for the first time. 

The predictions for neutron star masses as a function 
of the central density in Walecka, ZM and ZM3 models 
are shown in figure ^|. ZM model predicts a maximum 
mass of approximately 1.6Mq, in the limit of accept- 
ability for the mass of a pulsar. In particular, the ZM3 
model is very soft and predicts a very small maximum 
neutron star mass, ~ O.72Af . It may be surprising, 
at a first glance, that the maximum neutron star mass 
for the Walecka model with hyperons (2.77M©) exceeds 
the well known result (2.6M©) found in [Q for stars just 
composed of neutrons, since the addition of other par- 
ticles softs the equation of state, lowering the resulting 
star mass. This apparent contradiction can be explained 
by the extreme sensibility of this kind of theory on the 
specific choice of the values of the binding energy and 
saturation density. The authors cited above have used 
B = -15.75 MeV and p = 0.19 /m~ 3 ; with this choice 
we were able to reproduce their results. However, with 
our choice for these quantities, which is widely used in 
the recent literature, we get for the mass of a star com- 
posed only by neutrons the value of 3.05M©, that is, a 
difference of almost a half solar mass! Using the con- 
stants from |j (04 = 33.6 MeV), we get 2.33M© for the 
mass of a neutron star with the inclusion of hyperons and 
leptons. In this way, extrapolation for neutron star densi- 
ties from the fitting of B and po at saturation needs more 
precision on the choice of the values for these quantities. 

Figure [l3] shows the behaviour of the chemical poten- 
tials and field intensities in the Walecka, ZM and ZM3 
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models. These results should be compared to the corre- 
sponding one obtained with the BB model. We observe 
the same saturation of the electron chemical potential 
at ~ 200 M eV for the ZM and BB models with univer- 
sal coupling (see figure |l].d). ZM3 model presents higher- 
values for —g B go3, as compared to the other models; this 
behaviour is connected with the large value for g e /m e in 
this model (cf. tab le H). The known problem of nega- 
tive effective mass ]26|-|28| manifests itself dramatically 
in panel a for the Walecka model. We discuss this point 
in the next subsection. 

In figure |l4| the populations in a system consisting 
of hyperons, nucleons and leptons, for the Walecka and 
ZM models are shown. The poor description of neu- 
tron star phenomenology in the ZM3 approach lead us 
to omit its results for the baryonic and leptonic popu- 
lations. Walecka's baryonic distribution stabilizes after 
p ~ 1.0 fm and all species appear up till p ~ 0.7 fm 
which is approximately the density where \S\ exceeds M. 
The lepton populations never vanish in the ZM distri- 
bution and even at p ~ 1.2/m -3 baryonic species are 
emerging. Essentially, these differences are due to the 
strength of the scalar potential in these two models. 
Since a particle is created only when 



In both S and S-V cases, to a higher maximum neutron 
star mass it corresponds, in general, a less compressible 
matter (higher K) . To higher values of the nucleon effec- 
tive mass it corresponds weaker scalar potentials; since 
the latter is also directly related to the incompressibil- 
ity of nuclear matter, the maximum neutron star mass 
decreases with M*. In spite that the two cases repre- 
sent different descriptions of the neutron star problem, 
the figure shows that to a fixed value of the compression 
modulus it corresponds very close values for the neutron 
star mass; opposite to that, to different values of the star 
mass it correspond the same value for the nucleon ef- 
fective mass. This result indicates that the predictions 
of neutron star masses based mainly on the compression 
modulus are more model independent than those based 
on the nucleon effective mass. 

From figure 18 and table fv|, we see that for A < 0.5 
in the S-V case, we have obtained negative values for the 
nucleon effective mass, which corresponds to a density 
region for which the strength of the scalar condensate 
exceeds the free mass of the nucleon, \S\ > M. 



D. Negative effective mass 



qBfJ-n-qe,BfJ-e>gtB^0+gpBQ03l3B + (M B ~ g* trB (J ) 



(52) 



a large scalar field favors the early emergence of parti- 
cles. (One should compare theses results with the corre- 
sponding ones shown in section [II for the BB model.) In 
figure |l5| the predictions of the ZM model for the radial 
distribution of the different lepton and baryon species is 
presented. 

We analyze now our new class of models allowing vari- 
ations of the A parameter. Tables [V and [v| show results 
for the radius, redshift (z) and hyperon/baryon ratio for 
the maximum mass of a neutron star sequence, consid- 
ering hyperons, nucleons and leptons degrees of freedom, 
as a function of A; we present some nuclear matter prop- 
erties as well. In figure 16 we show the dependence of the 
maximum mass with this parameter. One can see that 
some models corresponding to the S-V case (including 
ZM3 model discussed above) predict very small neutron 
star masses, lower than the masses of all pulsars found 
until now. We observe again the similarity of the predic- 
tions associated to the ZM and ZM3 models compared to 
those with exponential couplings, EXP/S and EXP/S-V. 

The dependence of the maximum neutron star mass 
of a sequence with the nuclear matter bulk properties 
K and M* at saturation density is also presented in ta- 
bles [V and [v|. Figure |l^ exhibits the dependence of the 
maximum mass of a neutron star for a sequence with the 
compression modulus and nucleon effective mass at sat- 
uration density in S and S-V cases. Figure |l^ shows the 
maximum neutron star mass for a given value of the A 
parameter as a function of the scalar potential S in the 
center of the neutron star. 



The nucleon effective mass is a dynamical quantity and 
expresses the screening of the baryon masses by the scalar 
meson condensate. By analyzing in our approach the 
general expression for the baryon effective mass 



Ml = M B - 



g<j<j 



{l+g a a/{XM)Y 



(53) 



in the limit g a a — > oo, we see that only in case A > 1 
this quantity do not vanish or become negative. Some 
internal constraints in the theory can avoid this, as is the 
well known case of the Walecka model without hyperons 
M. However, as we add more and more baryonic species 
we open the possibility for the scalar potential \S\ to 
become greater than the free masses. Let us introduce 
the scalar field equation in ( |53] ) , for the Walecka case, to 
get a better understanding of the problem; we obtain 



M%=M B -J2 



k 2 dk 



-MtS 



M%: 




9o_ 
m„ 



k 2 dk 



k 2 



(54) 



M%< 



As we add more and more baryonic species, the negative 
term in the numerator of (|54|) becomes more important 
and we open the possibility for the scalar potential \S\ 
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to become greater than the free baryon masses. In fact, 
even if we had just considered nucleons, but taking into 
account the difference in the neutron and proton masses, 
that negative term in the numerator would appear. How- 
ever, since this difference is very tiny, a negative effective 
mass would emerge in practice only with the addition of 
hyperon degrees of freedom |^8j . 

We could interpret the vanishing of the effective mass 
as a signal of a transition to a quark-gluon plasma 
phase. However, we should remember that our La- 
grangian model does not comprehend these underlying 
degrees of freedom. Additionally, as pointed out in |30| , 
at such high densities and strong meson fields we have al- 
ready reached the critical density where the production 
of baryon-antibaryon pairs is favored. In fact, this behav- 
ior of the effective mass may indicate that the mean field 
approximation is reaching the limits of its applicability. 

Concerning the problem of negative effective mass, we 
want to clarify a point which may appear naive but has 
led to misunderstandings. The integrals appearing in the 
RMF formalism always present a term M* 2 and, in this 
way, are symmetric with respect to positive or negative 
values of the effective mass. For example, let us take the 
integral related to the scalar density: 

I(t,m) = f f* fc ; (55) 

JO V K + TCI 

rigorously, the expression for I{t, m) would involve the 
modulus of m since V m 2 = |m|: 

f* k 2 2, ft + y/F+m*} 

/ ==-t\/t z + m—-m In : — : . 

Jo Vfc 2 + m 2 2 v 2 \ \m\ J 

(56) 

Some authors have assumed ad hoc the modulus of M* in 
the above logarithm when the problem of negative values 
appears. However, we can see from expression (^6|) that 
it emerges naturally from the symmetry of the integrals. 
Of course, different results would arise if we are using 
\M*\ instead of M* in other expressions, e.g. the cubic 
term in the energy density expression for the BB model, 
since it can be rewritten as 

-bM(M-M*f. (57) 

Finally, we want to stress out that, in spite of the interest- 
ing issue about the physical interpretation of a negative 
effective mass, mathematically the RMF model continues 
to work even with M* < 0. 



VI. SUMMARY AND CONCLUSIONS 

We have done an analysis on the influence of nuclear 
matter properties on the structure of neutron stars using 
a new class of RMF models with parameterized couplings 



among mesons and baryons. These couplings allow to 
reproduce, trough a suitable choice of mathematical pa- 
rameters (A, /?, 7), results of the Walecka and derivative 
coupling models such as ZM and ZM3 models. As we 
have shown above, this new class of relativistic models 
permits some control on the intensity of the scalar, vector 
and, isovector mesons mean-field potentials. In particu- 
lar, we have scanned values of nuclear matter properties 
which correspond to intermediate region of predictions of 
the Walecka, ZM and ZM3 models. 

In this new class of models we considered two cases, the 
Scalar (case S) with variations of A and (3 = 7 = ( this 
case contains the results of the Walecka and ZM models) 
and the Scalar- Vector (case S-V) with variations of A 
and j3 = 7 = A (ZM3 model belongs to this category). 
Taking the limit A — > 00 for the S case, we have a model 
with exponential couplings (EXP/S model). Similarly, 
for A = /3 = 7->oowe have the asymptotic model 
EXP/S-V. 

An extended RMF which includes hyperon and lepton 
degrees of freedom was presented. For the sake of sim- 
plicity, we have assumed an universal coupling for the 
meson-baryon couplings. We have investigated the ef- 
fects of the non-linear coupling in the nucleon effective 
mass and compression modulus of nuclear matter and in 
neutron star properties. 

Having this in mind let us summarize the main findings 
of our work: 

• For the nuclear matter properties we have found 
values of parameters that reproduce reasonably the 
best range of phenomenological results, namely 0.6 < 
M*/M < 0.7 and 200 < K < 300 MeV. We have also 
found that to higher values of the nucleon effective mass 
it corresponds lower values of the compression modulus. 
The Scalar case leads to a narrow interval of values of 
A with different descriptions of nuclear matter proper- 
ties while the Scalar- Vector case has a broader range. To 
be more specific, case S gives reasonable results in the 
range A ~ 0.05 — 0.07 and case S-V better results for 
A ~ 0.16 — 0.4. We have also pointed out the strong sim- 
ilarity between the results of ZM-like models and those 
with exponential couplings. We have also obtained an 
analytical expression for the compression modulus of nu- 
clear matter as a function of the nucleon effective mass 
(Appendix §). 

• We have presented explicitly the neutron star phe- 
nomenology for the Walecka, ZM and ZM3 models. To 
the best of our knowledge, we display in this work, for 
the first time, the results of these two later models. The 
ZM model predicts a maximum mass of approximately 
1.6M Q for a neutron star, while the ZM3, being too soft, 
leads to ~ 0.72M Q as the limiting mass. We have ob- 
tained Walecka and ZM fcrmionic distributions. The cor- 
responding particle populations exhibit expressive differ- 
ences essentially due to the strength of the scalar poten- 
tial in these two models. 

• Sensibility of this approach to the specific choice of 
B and po was noticed with differences of the order of a 
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half solar mass. For the Walecka model, the neutron star 
mass with hyperons has varied from 2.77M Q to 2.33M Q 
with a very small change in those nuclear matter bulk 
properties. Thus, extrapolation for neutron star densities 
from the fitting of B and po at saturation needs more 
precision on the choice of the values for these quantities. 

• We analyzed our new class of models allowing the 
variation of the A parameter and showed results for the 
radius, redshift (z) and hyperon/baryon ratio for the 
maximum mass of a neutron star sequence, considering 
hyperon, nucleon and lepton degrees of freedom. In fig- 
ure |l^ we showed the dependence of the maximum mass 
with this parameter. We found that some models cor- 
responding to the S-V case (including the ZM3 model) 
predicted very small neutron star masses. We observed 
again the similarity between the predictions of neutron 
star phenomenology associated to the ZM and ZM3 mod- 
els compared to those with exponential couplings, EXP/S 
and EXP/S-V. 

• The dependence of the maximum neutron star mass 
of a sequence with the nuclear matter bulk properties 
K and M* at saturation density was also analyzed. In 
both S and S-V cases, to a higher maximum neutron 
star mass it corresponds, in general, a less compressible 
matter (higher K). We also concluded that the maximum 
neutron star mass decreases with the nucleon effective 
mass M*. 

• Finally, we found that for A < 0.5 in the S-V case, 
we have obtained negative values for the nucleon effective 
mass, which corresponds to a density region for which the 
strength of the scalar condensate exceeds the free mass 
of the nucleon, |5| > M. We have shown how the dif- 
ferences in the baryonic bare masses explain this result; 
as explained in the text, in practice the inclusion of hy- 
perons is responsible for this behaviour. These results 
indicate the existence of very intense scalar condensates 
in the interior region of neutron stars. Concerning math- 
ematics, the theory continues to work due to the symme- 
try of negative and positive values of M* in the Fermi 
integrals. 



We shall consider infinitesimal variations of the Gibbs 
potential (G = G(p,T, Ni)), with respect to the number 
of particles (Ni), at constant temperature T and pressure 
P 



dG = J^(— 



\dN 



dN, 



(A2) 



At chemical equilibrium, Gibbs energy obeys the condi- 
tion: 



E 



dG 
dN 



dN± 
dN 



0. 



(A3) 



with the ratio -r^ determined by the corresponding 
chemical reactions. Accordingly, if an element i suffers a 
variation i>i, the remaining elements will suffer a variation 
(vijvi)Vj to maintain the stoichiometry of the reaction. 
As a result, we have dNj/dN — v jl v i and we can write 
the condition of chemical equilibrium: 



^2 v iVi = 0, 



(A4) 



where the chemical potential of element i, fit, is defined 
as: 



dG 

m 



(A5) 



T,P,N i7 y 



In this way we observe that the chemical potentials obey 
the symbolic equation (Al), with the substitution of Aj 
by m. 

In general, if a chemical reaction respects given con- 
servation laws, the number of independent chemical po- 
tentials is equal to the number of these laws. 

In the following we consider two conservation laws, 
electric charge and baryon number. In this case we can 
express these laws, for a chemical reaction, as: 



N 



N 



^ Viq ei = and 2J v iQbi = 0, 



(A6) 
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APPENDIX A: CHEMICAL EQUILIBRIUM 

Chemical reactions, e.g. A + B ^ C + D, can be ex- 
pressed in general, as a symbolic linear combination of 
its components [511: 



(Al) 



For example, in the reaction n ^ p + + e + v we have 
for the coefficients v n = — v v = — v e = — v v = 1. 



where q e i and denote, respectively, the electric and 
baryon charges of element i. As we have, in this case, N 
variables and two equations, we are able to express only 
two coefficients Ui in terms of the remaining TV — 2, which 
will be independent: 



N 

v\q e \ + v 2 q C 2 = - y j Wei, 
i ^1,2 
N 

v\qb\ + Vi<lb2 = - 51 viqu- 
i ^1,2 



(A7) 



(A8) 



We consider as an example element 1, the neutron, and 
element 2, the electron. We then have qbi = —q e 2 = 1, 
<7h2 = <7ei = and the above equations become 
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N 
N 



(A9) 



Replacing (A£) in (A4) wc find: 



N N N 

^ ^Mi = X! (PnQbi)vi - ^2 {\l e q_ei)Vi. (A10) 

Since the i>\ are independent, the equality of this equation 
will be verified only if the coefficients are equal. From this 
expression then results 



Pi QbiPn Qeipe- 



(All) 



APPENDIX B: COMPRESSION MODULUS 



From the definition of the compression modulus, 



K = 



,d 2 {e/p) 



we have 



K 



, 2 / p d 2 P 
V \pdk 2 



dk 2 



(Bl) 



Q=0O 



P 



2p 1 dp 
P dp 



(B2) 



At saturation density, the pressure of nuclear matter is 



null. From expression (B2) we then have 



K = 3fc 



dp 



(B3) 



The chemical potential is a function of k and a. Thus, 
its derivative with respect to k is: 



dp dp da dp 
dk da dk dk 



(B4) 



From equation fl47|) for the chemical potential we have 
for the last term of expression (B4) 



9l 2fc 2 



k 



dp = m *2 

dk ml TT 2 ^ ^fc2 + M *2 



(B5) 



and for the derivative of the chemical potential with re- 
spect to the a field we get: 



dp 

da~ 



= 2m 



t dm *0 2fc 



M* 



8M* 



13 ml da 3tt 2 ^ ' k 2 + M *2 da 



(B6) 



To find the contribution of the da / dk term in expression 
(B4), we consider a profile function / = f(a, k), and the 
condition 



/(a, k) = m 2 a a + F{a)p s + G(a)m*„^p z = ; (B7) 



di 2 



'V 2 



from this expression 

da ~{df/dk) a 



with 



and 
da 



df 
dk 



dk (df/da) k ' 

T?t ^ d Ps . r,ni n * 9u d P 
~dk ^m?,^'dk ' 



(B8) 



(B9) 



dF(a) d Ps 

p s + F(a)— h m 



dG(a) 2 
— r nin — - — p 

da H ml da 

(BIO) 



da 

G(a)dg* p2 
ml da 



According to equations (|44|) and ([15|), the derivatives of 
the F(a) and G(a) functions are 



dF{a) 
da 



M 



\ + l g a a „*+2 



and 



dG{a) 2 



(3 + 1 1 ,2+2 



(Bll) 



(B12) 



To obtain an equation for the compression modulus of 
nuclear matter as a function of g a /m a and g^jm^, we 
have to consider the following ratios 

F(a) G { a) 1 8F(a) ^ ^d^a) (m3) 



9<r ' 9a9^ ' g 2 - 9cr 



9u>gl da 



Combining (|B5|), (|B6|), (gg) and (JBTOJ) with (g), we 
finally obtain: 



X = K, + K„ 4 KsA X KaB 



where 



K, 



K. 



2 



K 3A = 



K, 



-3K 



3k Q 



k 2 + M* 2 



(B14) 

(B15) 
(B16) 



fG(a)\ 2*g 



/c 2 + M* 2 v 9a 



13 V ffeff" 7 37T 2 

F(aY 



(B17) 
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K 3B = 



ml 



1 + 



F(a) 



2k£M* 



9a ' ^Jkl+M* 2 



2ml 










F(a) 2k*M* 




9- ^Ag+M* 2 



(B18) 



and 

K. 



= 1+1 4- 




(B19) 



[1] A. Hewish, S.J. Bell, J.D.H. Pikington, P.F. Scott, R.A. 

Collins, Nature 217, 709 (1968). 
[2] W. Baade and F. Zwicky, Phys. Rev. 45, 138 (1934). 
[3] S.L. Shapiro, S.A. Teukolsky, Black Holes, White Dwarfs 

and Neutron Stars (John Wiley, New York, 1983). 
[4] J.D. Walecka, Advances in Nuclear Physics 16 (Plenum 

Press, New York, 1986). 
[5] N.K. Glendenning, Astrophys. J. 293, 470 (1985). 
[6] J. Boguta and A.R. Bodmer, Nucl. Phys. A292, 413 

(1977). 

[7] J. Zimanyi, S. A. Moszkowski, Phys. Rev. C42, 1416 
(1990). 

[8] A. Delfino, C. T. Coelho, M. Malheiro, Phys. Rev. C51, 

2188, (1995); ibid Phys. Letts. B345, 361, (1995). 
[9] A. Delfino, M. Chiapparini, M. Malheiro, L. V. Belvede- 
re, and A.O. Gattone, Z. Phys. A355, 145, (1996). 

[10] M. Chiapparini, A. Delfino, M. Malheiro and A. Gattone, 
Z. Phys. A357, 47, (1997). 

[11] A. R. Taurines, Neutron Stars in relativistic mean field 
theories, dissertation (in Portuguese) (Universidade Fed- 
eral do Rio Grande do Sul, Porto Alegre, 1999) (unpub- 
lished) . 

[12] N.K. Glendenning, Compact Stars (Springer- Verlag, New 
York, 1997). 

[13] S.A. Moszkowski, Phys. Rev. D9, 1613, (1974). 

[14] B. K. Harrison, K.S. Thorne, M. Wakano, and J. Whee- 
ler, Gravitation Theory and Gravitational Collapse (Chi- 
cago University Press, Chicago, 1965). 

[15] JW. Negele and D. Vautherin, Nucl. Phys. , A207, 298, 
(1973). 



[16 
[17 

[18 

[19 

[20 
[21 

[22' 

[23; 
[24; 

[25 
[26; 
[27 
[28; 
[29 
[30 
[31 



R. C. Tolman, Phys. Rev. 55, 364 (1939). 

J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 

374 (1939). 

B.M. Waldhauser, J.A. Maruhn, H. Stocker, and W. 
Greiner, Phys. Rev. C38 1003 (1988). 
N. K. Glendenning, F. Weber, and S. A. Moszkowski, 
Phys. Rev C45, 844 (1992). 

M. Barranco et al, Phys. Rev. C44 178 (1991). 

A. Delfino, J. Dey, M. Dey and M. Malheiro, Phys. Letts. 

B363, 17 (1995); ibid Z. Phys. C51, 2188 (1999) 

Z. Qian, H. Song, and R. Su, Phys. Rev. C48, 154, 

(1993). 

M. Malheiro, A. Delfino, and C. T. Coelho, Phys. Rev. 
C58, 426 (1998). 

W. Koepf, M. M. Sharma, and P. Ring, Nucl. Phys. 
A533, 95, (1992). 

A. Bhattacharyya, S.K. Ghosh, Int. Journ. Mod. Phys. 
E7, 495, (1998). 

P. Levai, B. Lukacs, B. Waldhauser, and J. Zimanyi, 
Phys. Lett. B117, 5, (1986). 

R. Knorren, M. Prakash, and P.J. Ellis, Phys. Rev. C52, 
3470, (1995). 

J. Schaffner, and I.N. Mishustin, Phys. Rev. C53, 1416, 
(1996). 

R.J. Furnstahl, B.D. Serot, and H.-B. Tang, Nucl. Phys. 
A598, 539, (1996). 

I.N. Mishustin, L.M. Satarov, J. Schaffner, H. Stocker, 
and W. Greiner, J. Phys. G19, 1303, (1993). 
L. Landau, and E.M. Lifshitz Physique Statistique (Mir, 
Moscow, 1967). 



13 



TABLE I. Values of coupling constants, nucleon effective mass, compression modulus of nuclear 
matter and scalar and vector potentials at saturation density (ZM and ZM3 models). 



Model 


(ff CT /m CT ) 2 


(gjm^f 


(g P /m p ) 2 


M*/M 


K 


S 


V 




(fm 2 ) 


(fm 2 ) 


(fm 2 ) 




(MeV) 


(MeV) 


(MeV) 


ZM 


7.94 


2.84 


5.23 


0.85 


224 


-140 


84 


ZM3 


19.57 


13.45 


9.06 


0.71 


159 


-267 


204 



TABLE II. Nucleon effective mass and compression mod- 
ulus of nuclear matter for different type of couplings between 
the scalar mesons and the nucleon fields |24|. 



m*(a) 


M*/M 


K(MeV) 


1 1 SctO- 
x M 


0.55 


545 


2 l-tanh(^cT) 


0.71 


410 


3 exp(-Zff) 


0.80 


265 


4 (l + ^ff)- 1 


0.85 


233 


TABLE III. Values of A, /3 and 7 for different QHD models. 


Model A 





7 


Walecka 








ZM 1 








ZM3 1 


1 


1 
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TABLE IV. Stellar properties for the S case: e c - central density;M* star mass; 71* - star radius; S - scalar potential in 
the star center; z - redshift; Y/A is the hyperon/baryon ratio and Nbt is the total baryonic number. All this quantities are 
evaluated for the neutron star with the maximum mass in the sequence. Besides we have the compression modulus K, the 
nucleon effective mass M* /M and the relativistic coefficient R. 



A 


log(e c ) 


M* 


R* 


S 


z 


Y/A 


Nbt 


K 


M*/M 


R 




I 3 

g/cm A 


(M & ) 


(km) 


(MeV) 






(x 10 ) 


(MeV) 









15.18 


2.77 


13.17 


936 


0.623 


0.27 


0.40 


566 


0.537 


0.931 


0.03 


15.24 


2.56 


12.39 


934 


0.597 


0.30 


0.36 


396 


0.598 


0.943 


0.05 


15.31 


2.35 


11.63 


929 


0.574 


0.32 


0.33 


310 


0.650 


0.951 


0.07 


15.38 


2.17 


10.89 


923 


0.554 


0.34 


0.30 


258 


0.694 


0.957 


0.09 


15.43 


2.02 


10.38 


910 


0.533 


0.35 


0.28 


235 


0.725 


0.960 


0.11 


15.47 


1.91 


9.98 


896 


0.515 


0.36 


0.26 


223 


0.749 


0.962 


0.13 


15.49 


1.83 


9.75 


877 


0.495 


0.35 


0.25 


217 


0.766 


0.964 


0.15 


15.52 


1.77 


9.59 


857 


0.479 


0.35 


0.24 


216 


0.779 


0.965 


0.17 


15.54 


1.72 


9.45 


838 


0.467 


0.35 


0.23 


213 


0.789 


0.966 


0.20 


15.53 


1.68 


9.48 


807 


0.446 


0.33 


0.23 


212 


0.798 


0.967 


0.25 


15.53 


1.61 


9.49 


732 


0.416 


0.30 


0.22 


212 


0.814 


0.968 


0.30 


15.52 


1.59 


9.61 


669 


0.399 


0.27 


0.21 


214 


0.822 


0.969 


0.35 


15.51 


1.58 


9.69 


618 


0.389 


0.26 


0.21 


216 


0.828 


0.969 


0.40 


15.51 


1.58 


9.73 


580 


0.385 


0.25 


0.21 


218 


0.833 


0.969 


0.60 


15.49 


1.58 


9.86 


480 


0.377 


0.22 


0.21 


223 


0.843 


0.970 


1.00 


15.47 


1.59 


9.98 


401 


0.372 


0.20 


0.21 


224 


0.850 


0.970 


1.50 


15.47 


1.59 


9.98 


366 


0.373 


0.20 


0.21 


226 


0.854 


0.971 


OG 


15.47 


1.59 


10.00 


350 


0.373 


0.20 


0.21 


228 


0.856 


0.971 






TABLE V. 


Stellar properties for the S-V case. 


Same correspondences 


as in table 


IV 




A 


log(e c ) 


M* 


R* 


S 


z 


Y/A 


Nbt 


K 


M*/M 


R 




i 3 

g/cm 


(Mo) 


(km) 


(MeV) 






(xlO 58 ) 


(MeV) 









15.18 


2.77 


13.17 


936 


0.623 


0.27 


0.40 


566 


0.537 


0.931 


0.03 


15.20 


2.70 


12.93 


944 


0.615 


0.28 


0.39 


510 


0.545 


0.933 


0.05 


15.22 


2.63 


12.64 


954 


0.610 


0.29 


0.38 


458 


0.554 


0.935 


0.07 


15.24 


2.56 


12.39 


960 


0.602 


0.30 


0.37 


417 


0.561 


0.936 


0.09 


15.27 


2.50 


12.12 


969 


0.598 


0.31 


0.36 


387 


0.567 


0.938 


0.11 


15.28 


2.43 


11.89 


973 


0.588 


0.32 


0.34 


358 


0.574 


0.939 


0.13 


15.30 


2.37 


11.68 


977 


0.579 


0.33 


0.33 


339 


0.579 


0.940 


0.15 


15.33 


2.30 


11.38 


985 


0.574 


0.34 


0.32 


311 


0.587 


0.941 


0.17 


15.35 


2.24 


11.16 


986 


0.563 


0.34 


0.31 


293 


0.594 


0.942 


0.20 


15.38 


2.17 


10.88 


993 


0.559 


0.35 


0.30 


276 


0.600 


0.943 


0.30 


15.51 


1.83 


9.58 


1011 


0.516 


0.39 


0.25 


218 


0.630 


0.949 


0.35 


15.56 


1.70 


9.09 


1012 


0.491 


0.41 


0.23 


205 


0.640 


0.950 


0.40 


15.62 


1.57 


8.60 


1014 


0.470 


0.42 


0.21 


195 


0.649 


0.951 


0.60 


15.62 


1.07 


8.08 


891 


0.282 


0.35 


0.14 


169 


0.682 


0.955 


1.00 


15.31 


0.72 


9.76 


577 


0.128 


0.10 


0.09 


159 


0.710 


0.959 


1.50 


15.18 


0.67 


10.21 


468 


0.113 


0.04 


0.08 


156 


0.728 


0.961 


CO 


15.14 


0.66 


10.31 


431 


0.110 


0.03 


0.08 


155 


0.738 


0.961 
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FIG. 2. Panel a: distribution of the energy density throughout a star composed of nucleon s an d leptons. Panel b: Neutron 
star mass as a function of the central density for matter with (curve I corresponds to x = \/2/3 and curve II to \ — 1) an d 
without (curve III) hyperons. Panel c: mass- radius relation (same labels as in panel b). 
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6). 



FIG. 5. Evolution of the coupling constants g% lU /m atU> with increasing density for ZM model [panel a) and ZM3 model [panel 
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FIG. 6. Dependence of the coupling constants with the A parameter for the scalar (full lines) and scalar-vector (dotted lines) 
models. 
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FIG. 7. Plane of the scalar S and vector V potentials. The potentials in the Scalar and Scalar- Vector exponential coupling 
models are indicated." 
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FIG. 8. Compression modulus of nuclear matter K as a function of the A parameter. 
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FIG. 9. The ratio M* /M, at saturation density, as a function of the A parameter. 
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FIG. 10. Comparison of the compression modulus of nuclear matter K with the ratio M* /M for the S and S-V cases (the 
values for the exponential models are indicated). The box shows the range of accepted values. 
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FIG. 17. Dependence of the maximum neutron star mass of a sequence with the compression modulus (left) and nucleon 
effective mass at saturation (right). Solid line corresponds to S case and dashed line to S — V case. 
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FIG. 18. Maximum neutron star mass as a function of the scalar potential at the star center. Solid line corresponds to S 
case and dashed line to S — V case. The long dashed line corresponds to S = 939 MeV. 
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